Valence bond glass — A unified theory of electronic disorder and pseudogap phenomena 

in high T^ cuprate superconductors 
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We show that the low-energy fluctuations of the valence bond in underdoped high-T^ cuprates, originating 
from quantum fluctuations of the superexchange interaction, are pinned by the electronic disorder due to off- 
stoichiometric dopants, leading to a valence bond glass (VBG) pseudogap phase. The antinodal Fermi surface 
sections are gapped out, giving rise to a normal state Fermi arc whose length shrinks with underdoping. Below 
Tc, the superexchange interaction induces a J-wave superconducting gap that coexists with the VBG pseudogap. 
The evolution of the local and momentum- space spectroscopy with doping and temperature captures the salient 
properties of the pseudogap phenomena and the electronic disorder. The unified theory elucidates the important 
interplay between strong correlation and the intrinsic electronic disorder in doped transition metal oxides. 



The least understood electronic state in the high-T^ cuprates 
lies in the pseudogap phase that straddles the AF Mott insu- 
lator and the J-wave superconductor. The origin of the pseu- 
dogap phase is a subject of intensive debate [1]. This unex- 
pected phase is defined by a normal state pseudogap in the 
single-particle excitation spectrum below a characteristic tem- 
perature r* 1I2II3I]. One of the most intriguing properties is the 
momentum space anisotropy: the Fermi surface (FS) near the 
J- wave antinodes is gapped out, leaving a Fermi arc of gapless 
excitations near the nodes. This nodal- antinodal dichotomy is 
further enlightened by recent ARPES ||1[10] and Raman ^ 
experiments that find two distinct spectral gaps in terms of 
the temperature and doping dependence: a large pseudogap 
near the antinodes that tracks T* and a smaller superconduct- 
ing (SC) gap along the arc that tracks Tc. Recent STM mea- 
surements observe a SCgap coexisting spatially with a large 
pseudogap below T^ llsl ISD • 

The Fermi arc raises the possibility for an unconventional 
electronic state. In the absence of disorder, the FS is the trajec- 
tory of the poles in the single-particle Green's function (dis- 
continuity in the momentum distribution) and involves con- 
tinuous contours in momentum space. Thus the Fermi arc 
must evolve into either a Fermi point or a FS pocket at low 
temperatures. This difference highlights two different (one- 
gap and two-gap) proposals for the origin of the pseudogap. 
In the one-gap scenario, the pseudogap is a J-wave pairing 
gap without SCphase coherence that evolves into a single SC 
gap below Tc ||ll liol luJ] . The two-gap scenario attributes the 
pseudogap as due to competing order from an energetically 
favorable state unrelated to pairing. An important implication 
of a generic two-gap scenario is the coexistence of a J- wave 
SC gap over the Fermi arc and a low-temperature pseudogap 
in the antinodal region below Tc [12]. The observation of two 
distinct gaps by ARPES, Raman, and STM experiments pro- 
vides strong support for the two-gap scenario. 

A consistent theory for the pseudogap has been challeng- 
ing. Most of the proposals for a competing state invoke var- 
ious density waves and flux/orbital current order ill llOB that 
break translation symmetry by one lattice constant and result 
in superstructures. The FS is truncated in the antinodal region 
by band-folding with respect to a unique ordering wavevector 



q, leading to FS pockets and particle-hole asymmetric pseu- 
dogap density of states lll3ll . However, despite many years 
of effort and the much improved measurement resolution, no 
signatures of disguised FS pockets or folded bands have been 
detected in zero external fields. 

A related phenomenon in the cuprates is the electronic 
disorder. There has been mounting STM evidence for 
nanoscale DOS gap disorder Q [ll H [HI Q and, 
under a variety of conditions where superconductivity is 
weakened, short-range ordered checkerboard DOS modula- 
tions lllSl [18, 19, 20] which are manifestations of a bond- 
centered electronic glass 112 111 . A natural cause for the elec- 
tronic disorder is the out of plane ionic dopants, intersti- 
tial in Bi2Sr2CaCu208+x. substitutional in La2-xSrxCu04 and 
Ca2-xNaxCu02Cl2, and in combination with chemical substi- 
tutions in Bi2Ln2-2Bi^Cu06+x (Ln-Bi2201). In addition to 
inducing structural distortions, the screening of the dopant 
electrostatic potential is highly nonlinear in doped Mott in- 
sulators due to strong Coulomb repulsion, leading to inhomo- 
geneous electronic states with spatial variations in the local 
doping concentration ||22l |23|] . 

We present here a unified theory for the pseudogap phe- 
nomena and the electronic disorder. The key point is the 
interplay between strong correlation induced valence bond 
fluctuations and the dopant induced disorder. In the parent 
compounds, the superexchange interaction between the Cu 
spins, described by the nearest neighbor Heisenberg model 
Hj = JY.{ij) ^i ' Sjs causes AF order in the ground state. Due 
to strong quantum fluctuations of the spin- 1/2 moment, the 
nonmagnetic valence bond states are close in energy to the AF 
state. Writing S^ • Sy = \c]^Ci(yic\^,Cj(j + const., the valence 
bond can be formed via spin-singlet pairing A^^ = {c^C j \ — 
C;|Cjl) as in the resonance valence bond (RVB) theory |[24ll . 
and the paramagnetic Xij = Lcr i^^Ja^jo) as envisioned by Paul- 
ing in chemical bonding. Since charge fluctuations are com- 
pletely suppressed, the two descriptions are equivalent and the 
valence bond states are highly degenerate owing to the SU(2) 
symmetry OSll . Besides valence bond liquid states, there are 
also symmetry breaking valence bond crystal states that are 
gapped but competitive in energy 111 |26|] . Doping the Mott 
insulator breaks the SU(2) symmetry and makes the valence 



bond in the particle-particle and particle-hole channels differ- 
ent. The basic question is which fluctuating valence bond state 
is selected when a sufficient amount of holes destroys the AF 
long-range order. In the short-range RVB theory, the spin- 
singlet pairs are mobilized by the doped holes and condense 
into a J- wave SC state [27]. There is, however, a natural com- 
peting order driven by the same superexchange interaction, 
i.e. the paramagnetic valence bond Xij- 

We show that the J-wave component of the real part of 
Xij represents the most important low-energy fluctuations and 
that dopant induced disorder pins such J- wave charge density 
wave fluctuations to an electronic valence bond glass (VBG) 
phase exhibiting the observed pseudogap phenomena. This 
is similar in spirit to disorder induced glassy phases or the ne- 
matic liquid crystal of stripes [28, 29]. A ubiquitous feature of 
the VBG is the emergence of a pseudogap near the antinodes 
where the valence bond fluctuations are large due to the flat- 
ness of the band, giving rise to a genuine normal state Fermi 
arc. Based on microscopic calculations of the extended t-J 
model using spatially unrestricted Gutzwiller approximation, 
we show that the VBG captures the salient properties of the 
observed pseudogap phenomena and electronic disorder, es- 
pecially in underdoped and chemically substituted bilayer and 
single-layer Bi-based cuprates where the pseudogap and the 
SC gap are wefl separated iH B 0, BS] . 

The Hamiltonian of the extended t-J model is given by 

^ = -"LUjPiclcjaPj^J Y. (Si • Sj - \hihj)^Y.^ihi. (1) 
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The electrons hop between near neighbors via ttj. Repeated 
spin indices are summed, fii = cJ^Cia, ^i = ^<^Ja^a&<^i&- The 
projection operator Pi removes double occupation on site-/. 
The last term in Eq. ([T]) is the electrostatic potential, V{i) = 



^JT^i \n-rj\ W 
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Coulomb interaction of strength Vc between the in-plane elec- 
trons provides the important screening of the ionic potential 
of strength Vj from .^number of off-plane dopants at a set 
back distance J^ E |23l] . We use / = 120meV and up to fifth 
nearest neighbor hoppings t = (360,-120,29,24, — 24)meV 
relevant for the band structure ll30l] and set V^ = Vj = 0.5eV 
and J^ = 1 in units of the lattice constant |22]. To account 
for their Coulomb repulsion, the ionized dopant configura- 
tions are generated randomly with a hard-core of one to three 
lattice spacings. 

To account for both strong correlation and disorder, we use 
the spatially unrestricted Gutzwiller projected wave function 
1^) = ni\yr(l ~%%)l^o). where ^o is a Slater determi- 
nant state and yi is a local fugacity that keeps the density un- 
changed before and after projection. The projection is imple- 
mented using the Gutzwiller approximation [Sll [321 • The t-J 
Hamiltonian is replaced by one in the unprojected space with 
renormalized hopping and exchange that capture the basic 
Mott physics: Uj -^ g\:tij and / 
state, g\j = ^AxiXj/(l^Xi){l- 



gjjj. In an inhomogeneous 
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Figure 1: Normal state pseudogap and Fermi arc at x = 0.125. (a) 
Distribution of J- wave valence bond showing glassy order, (b) Aver- 
aged DOS showing VBG pseudogap. (c) and (d) Fermi arc. 



depend on the local doping Xi llj, |23|] . Decoupling the ex- 
change interaction in terms of the valence bond Xij and singlet 
pairing A^y, we obtain a Gutzwiller renormalized Hamiltonian, 



^GA = - X g\jtijc]^Cjcy + X( V,- + h)c]^Cicy - ^ XiUi 

- ^-^ L 4 (Xv^^iCT +h.c. - \Xijf) 



(iJ) 



- -7jY.stj{^ij£aa'CiaCja'^^^C.-\Aij\^), (2) 

^ {iJ) 

where gf- = gf- = gj- and Xi originates from the local fugacity. 
We minimize the ground state energy of Eq. ([S]) through self- 
consistently determined {xi^Xi^Xij^^ij} on 24 x 24 systems 
for different dopant configurations. 

Normal state VBG pseudogap phase. In the normal state 
above Tc or the zero temperature phase with Aij = 0, we find 
that the valence bond Xij is real and fluctuates due to the disor- 
der potential in the doping range studied. The dominant fluc- 
tuation is in the J- wave channel, Xd{i) = \ Hj^ijXij^ with the 
form factor dij = ±1 for the four bonds emanating from site /. 
We first focus on the normal state by setting Aij = 0. The his- 
togram of Xd at average doping x = 0. 125 is shown in Fig. la. 
It follows a Gaussian distribution with zero mean. The root- 
mean-squared fluctuation represents a nonzero "glassy" or- 



der parameter 5^ = Jl^iXJi^/^s = 0.28 for the VBG phase. 
We note that glassy dynamics of valence bond was studied in 
the weak-coupling metallic phase of the Hubbard-Heisenberg 
model at half-filling 1330 . The most succinct feature of the 
VBG is the emergence of the pseudogap and the Fermi arc. 
The averaged DOS in Fig. lb shows a remarkable V-shaped 
pseudogap, A^^ ~ 32meV, approximately symmetrically dis- 
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Figure 2: Doping dependence of the Fermi arc and valence bond 
glassy order (a) and averaged DOS and pseudogap (b). (c) Tempera- 
ture evolution of VBG pseudogap at x = 0. 125. 
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Figure 3: Angule dependence of symmetrized spectral function (a) 
and pseudogap (b). (c) Doping dependence of pseudogap from DOS 
(circles) and J- wave valence bond order parameter (squares), (d) 
\Xd{Q)\ along ^;c-direction with peaks at q* marked in Fig. Id. 



tributed around the Fermi level due to the J- wave nature of the 
VBG. The calculated spectral intensity at the Fermi energy in 
Fig. Ic reveals a FS truncated near the antinodes by the pseu- 
dogap and a Fermi arc around the nodes. The Fermi arc tracks 
the underlying FS and remains prominent upon lowering the 
intensity scale in Fig. Id without signs of band folding. 

The doping and temperature dependence of the VBG pseu- 
dogap is shown in Fig. 2. With increasing doping, the dis- 
tribution of the valence bond in Fig. 2a sharpens around 
zero, the VBG order parameter 5^ reduces, and the Fermi 
arc length increases in agreement with ARPES experiments. 
Note that while the density of off-plane dopants increases 
with doping, 5^ decreases due to improved screening by more 
mobile carriers that leads to weaker fluctuations of the po- 
tential Vi. The doping evolution of the averaged DOS is 
shown in Fig. 2b. The pseudogap becomes smaller and shal- 
lower with increasing doping and becomes undiscernible be- 
yond X = 0.22 on the shoulder of the van Hove peak. The 
extracted pseudogap has a doping dependence that follows 

Apg - (//8) ^LiiXjglXijdijY/Ns. as shown in Fig. 3c. The 
temperature evolution of the pseudogap obtained by minimiz- 
ing the free energy of Eq. ^ is shown in Fig. 2c at x = 0. 125. 
The pseudogap onset temperature T* :^ 16meV is clearly seen 
to be determined by the thermal filling of the VBG pseudogap. 

The momentum-dependence of the symmetrized spectral 
function A(/:, (o) is plotted in Fig. 3a at Fermi level for x = 
0.14. It shows gapless quasiparticle excitations along the arc 
and the opening of the pseudogap at the arc tip that increases 
toward the antinode. The angular dependence of the pseudo- 
gap extracted this way is shown in Fig. 3b. Away from the 
zero-gap Fermi arc regime, the pseudogap follows the J- wave 
form (dashed line) consistent with the J- wave nature of the 



VBG. The line-shape in Fig. 3a shows that the pseudogap near 
the antinodes is a soft gap. This important prediction of the 
VBG is consistent with the original ARPES experiments \MM 
and differs from the naive picture of a hard gap that depletes 
all states near the antinode below the pseudogap energy scale. 
These in-gap states contribute to the spectral weight at low 
energy and have important consequences in the SC state. 

The pseudogap in the VBG theory is not sensitive to the 
precise topology of the FS and does not require the antin- 
odal FS sections to be parallel which can induce long-range 
incommensurate density wave order lll3n . We obtain simi- 
lar results using a different set of hopping parameters t = 
(480, -160,50,50, -50)meV used in Ref.lQ that does not 
favor the antinodal nesting condition. A band dispersion 
showing the observed flatness near the antinode suffices. The 
propensity toward the VBG is due to an enhanced static sus- 
ceptibility broadly peaked around the q" shown in Fig. 1 that 
enables electronic disorder to pin the VBG with a distribution 
of Xdi^) peaked around q". In Fig. 3d, we show the Fourier 
power spectrum \Xd{Q)\ plotted along the ^^-direction. Both 
the peak and the incommensurate q" decrease with increasing 
doping. At X = 0.14, ^* c^ (±0.4;r,0), indicative of glassy 
~ 5 X 5 checkerboard patterns for the valence bond. 

Superconducting phase. The VBG can coexist with an in- 
homogeneous J-wave superconductor; both due to the su- 
perexchange interaction. Although the Gutz wilier factors 
gfj = gfj^ charge fluctuations at finite doping, and in particu- 
lar, the pair-breaking induced by inter- site Coulomb repulsion 
will weaken the singlet pairing channel ll34ll . Moreover, the 
electron-phonon interactions have been shown to promote the 
J- wave charge density wave OSl l36l] . To incorporate these 
effects into the renormalized mean field theory in Eq. ([2l), we 



use gf- = g/. and gf 
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Figure 4: SC state 3.tx = 0.14. (a) and (b) Averaged DOS showing 
SC gap and VBG pseudogap. (c) Temperature dependence of A^c and 
Apg. (d) LDOS along two line cuts, (e) Angular dependence of SC 
gap and pseudogap. (f) Age and A^^ as a function of doping. 



the two energy gap scales in the underdoped regime. The av- 
eraged DOS at r = is shown in Fig. 4a for x = 0.14. It 
displays two gaps: a smaller SC gap A^^ — lOmeV and a large 
pseudogap A^^ c^ 22meV inherited from the normal state, in 
agreement with STM experiments on La-Bi2201 [8, 9]. The 
temperature evolution in Fig. 4b shows that the SC gap and 
coherence peaks disappears above Tc ~ 4meV as the system 
enters the pseudogap phase. Interestingly, as T is increased 
toward T^ the pairing gap and the pseudogap show opposite 
temperature dependence (Fig. 4c), a typical feature of coex- 
isting but competing order. 

In Figs. 4d, the local DOS (LDOS) is shown along two line 
cuts on a typical sample at x = 0.14 . The evolution of the 
line-shape agrees with the STM conductance spectra observed 
in La-Bi2201 ||8l |9|], exhibiting two coexisting low energy 
gaps that are spatially inhomogeneous. The momentum de- 
pendence of the two gaps, calculated from the single-particle 
spectral function A (^, ft)), is plotted along the underlying FS in 
Fig. 4e. Remarkably, the J- wave pairing gap extends beyond 
the Fermi arc into the antinodal regime. The coexistence of 
the two gaps off the Fermi arc in momentum space is an im- 
portant prediction of the present theory. Physically, this is 
a consequence of the softness of the normal state pseudogap 



discussed above which allows pairing of the antinodal states 
inside the pseudogap. The ground state below Tc is thus a co- 
herent mixture of J- wave VBG and SC pairs; there are traits 
of the glassy valence bond order in the pairing gap near the 
antinode and vice versa. Although incoherent background and 
inelastic life-time broadening tend to mask the coherent peak 
and the pairing gap near the antinode in the spectral function, 
recent high resolution ARPES experiments on La-Bi2201 in- 
deed observe two gaps near the antinodes 1191 [3711. 

In summary, we presented a VBG theory for the essential 
features of the normal state pseudogap and the two-gap phe- 
nomena in the SC state. In Fig. 4f, we construct a theoretical 
phase diagram using the doping dependence of the J-wave 
pairing gap and the VBG pseudogap. It captures the basic 
topology of the global phase diagram of the high-T^ cuprates. 
Although the present theory does not include a precursor pair- 
ing induced pseudogap above T^ it does not rule out such a 
possibility due to fluctuations beyond the Gutzwiller theory. 
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